Quantum theory of spectral diffusion induced electron spin decoherence 
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A quantum cluster expansion method is developed for the problem of localized electron spin 
decoherence due to dipolar fluctuations of lattice nuclear spins. At the lowest order it provides 
a microscopic explanation for the Lorentzian diffusion of Hahn echoes without resorting to any 
phenomenological Markovian assumption. Our numerical results show remarkable agreement with 
recent electron spin echo experiments in phosphorus doped silicon. 
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It was realized a long time ago that spectral diffusion 
due to the dipolar fluctuations of nuclear spins often dom- 
inates the coherence decay in electron spin echo exper- 
iments P, 01 ■ The recent advent of spin-based quan- 
tum computation in semiconductor nanostructures re- 
vived the interest in spectral diffusion, which is expected 
to be the dominant channel for low-temperature spin de- 
coherence in several spin-based quantum computer ar- 
chitectures 0. In spectral diffusion, the electron spin 
Zeeman frequency diffuses in time due to the noise pro- 
duced by the nuclear spin bath. Dipolar fluctuations in 
the nuclear spins give rise to a temporally random effec- 
tive magnetic field at the localized electron spin, lead- 
ing to irreversible decoherence (i.e. a T2-process). All 
available theories to date are based on classical stochas- 
tic modeling of the nuclear field, a Markovian theoretical 
framework which is inevitably phenomenological since it 
requires an arbitrary choice for the spectrum of nuclear 
fluctuations. Such a classical Markovian modeling is ar- 
guably incompatible with the strict requirements of spin 
coherence and control in a quantum information device. 
In addition, recent rapid experimental progress in single 
spin measurements 4], which in the near future promise 
sensitive measurements of quantum effects in spin reso- 
nance, also warrant a quantum theory of spectral diffu- 
sion. Here we propose a quantum theoretical framework 
for spectral diffusion which is non-stochastic and fully mi- 
croscopic. In addition, our theory produces an accurate 
quantitative prediction for the initial decoherence, which 
is the most important regime for quantum computation. 
To the best of our knowledge, ours is the first quantum 
theory for electron spin spectral diffusion. 

Spectral diffusion is not a limiting decoherence pro- 
cess for silicon or germanium based quantum computer 
proposals because these can, in principle, be fabricated 
free of nuclear spins using isotopic purification. Unfortu- 
nately this is not true for the important class of materials 
based on III-V compounds, where spectral diffusion has 
been shown to play a major role There is as yet no 

experimental measurement of localized spin decoherence 



(echo decay) in III-V materials, but such experimental 
results are anticipated in the near future. 

Our theory reveals that the inclusion of quantum cor- 
rections to nuclear spin fluctuation increases the degree of 
decoherence, as is best evidenced from our explanation 
of the existing factor of three discrepancy between the 
Markovian stochastic theory |5j| and experimental data 
0, 0, of spin echo decay in phosphorus doped silicon. 
Our method allows a fully microscopic explanation for 
the observed time dependence of Hahn echo decay due 
to a nuclear spin environment. It was pointed out a long 
time ago 00 ^^^^ the observed time dependence of these 
echoes are well fitted to the expression exp (— r^) (here r 
is half of the time lag between the initial signal and an 
echo), a behavior which can be derived phenomenolog- 
ically by assuming Lorentzian Brownian motion for the 
electron spin Zeeman frequency 00. In our method 
this behavior arises naturally from the collective quan- 
tum evolution of the dipolar coupled nuclei, without any 
phenomenological assumption on the dynamics of the en- 
vironment responsible for decoherence. A proper descrip- 
tion of coupled spin dynamics is rather difficult due to 
the absence of Wick's theorem for spin degrees of free- 
dom. In that regard, variations of our method may prove 
rather useful, since environmental spin baths are ubiqui- 
tous in any device exploiting the coherent properties of 
quantum spin systems. 

The free evolution Hamiltonian for the spectral diffu- 
sion problem is given by 0| 
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Here S denotes the electron spin operator which couples 
to the nuclear spin J„ located at the lattice site i2„ . The 
nuclear spins are coupled to the electron through the hy- 
perfine constant A„ . We have truncated Eq. Q since the 
non-secular hyperfine coupling can be neglected at mod- 
erate magnetic fields {B > 0.1 Tesla for the Si:P case) . 
This interaction leads to interesting efl^ects at B = [lOj , 
but at the moderate magnetic fields required for spin 
resonance measurements it only contributes a small vis- 
ibility decay Each nuclear spin is coupled to all 
others via the dipolar interaction Eq. ((S} , which is again 
truncated in the range of moderate B fields (For further 
details we refer to Ref. 0). The Hahn echo experiment 
consists in preparing the electron spin in the initial state 
IVe) — (IT) + * U))/V25 E^'^d then allowing free evolu- 
tion for time t. A 7r-pulsc (here described by the Pauli 
operator a^^e) is then applied to the electron spin, and 
after free evolution for one more interval t an echo is 
observed, which provides a direct measurement of single 
spin coherence (i.e. T2 or Tm in the usual notation). 

We will now derive an exact expression for the Hahn 
echo decay due to Eq. The density matrix (for elec- 
tron and nuclear spins) describing Hahn echo is given by 



p(t) = U{t)poU\t) 
with the evolution operator 
C/(r) 



(6) 



(7) 



Here po is taken to be a thermal state for the nuclear 
spins. 



P0 = ^|ye)(2/e|®e-«"/'=-^, 
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where Tin — 'Hzn +'Hb and M is its partition function 
(M « 2^ for T > nK 0, where N is the number of 
nuclear spins). The spin echo envelope is then given by 



VEiT)=2\TT{{S^+iSy)p{T)}\. 



(9) 



An explicit expression for Eq. ||5J) can be obtained by not- 
ing that the electron and nuclear spin Zeeman energies 
commute with the total Hamiltonian, and f7x,e anticom- 
mutes with Sz. After a few manipulations we get 
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ve{t) 
where 

?7±(r) =e-^^±^ (11) 
are evolution operators under the effective Hamiltonians 



(12) 



which describe dipolar evolution under the effect of an 
electron spin up (H+) or down {Ti-). The trace in 
Eq. H10|) is taken over nuclear spin states only. 

In the high temperature limit (kgT ^ liB) we can 
expand Eq. H10|) in powers of r to get 
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Defining the parameter 



Ah n.m. 
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we obtain the first five coefficients D21 as a power series 
of Cnm and bnm- For example, the first two coefficients 
become explicitly D2 = 0, D4 = ^Y.n<m. ^nmKm- Trun- 
cating Eq. (|13|l gives physical results only for extremely 
short T unless most nuclear pairs satisfy the condition 
Cnm 1- Nevertheless most physical problems are char- 
acterized by several c„„i ^ 1, making evident the need 
for an alternative expansion. Thus the t expansion, while 
being formally exact, is not practical for coherence cal- 
culation except for extremely short r. 

In the S> 1 regime, non-degenerate perturbation 
theory is applicable to Eq. (|12|l . We introduce a book- 
keeping parameter A such that ±.Ti± — Ho ± XH.' . Here 
the unperturbed Hamiltonian Hq — i Aninz is diag- 
onal in the nuclear spin z-basis, while TL' — \'Hb is the 
dipolar interaction rescaled to have the same magnitude 
as Tio- A convenient choice for an order of magnitude 
estimate of A ~ l/|c„„i| is to use the minimum possi- 
ble value of bnm /\ An — Am\ between nearest neighbors: 
A ^ Max{bnm)/Max{An) ~ 10"^ for Si:P. Below we in- 
troduce a cluster expansion that can be related to powers 
of A in this perturbation approach when Cnm 2> 1. 

Let V he a subset of the nuclei in the problem. Let 
vv{t) be the solution of ve{t) [Eq. (|10|l ] when only in- 
cluding the nuclei in V. We recursively define the contri- 
bution from subset V as v-p (t) minus contributions from 
any proper subset of I?, 



E -sir)- 
scv 
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For the empty set, we define f0(r) vijj^r) = 1. 
Consider a subset contribution written in the form of 
the non-degenerate perturbation expansion. Assuming 
Max(6„m)r <C 1, we can show by the specific properties 
of TC' that a cluster of size k is composed of terms that 
are 0{X'') or higher. In other words, we can write the 
following expansion. 



VEiT) = i + j2 E ^pM + o(a^''+') 

k=2 \V\=k 



(16) 



where the second summation is over all possible nuclear 
subsets of size k (containing k distinct nuclear sites) . We 
note that a subset of size 1 gives no contribution. 
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The 0{X''"~^^) error in Eq. H16|l is misleading because 
the number of terms of a given order of A may be large 
compared to A. The nature of this problem, as well as a 
solution, becomes apparent when we use a nearest neigh- 
bor approximation. With this approximation, we ignore 
the interaction between distant nuclei (i.e. pairs of nuclei 
for which bnm is below some threshold) and divide our 
nuclear subsets into connected "clusters." A subset con- 
tribution is then the product of its cluster contributions. 

Consider all possible contributing subsets of size k. We 
can categorize these subsets by the number of clusters 
they contain. To estimate the number of subsets of size 
k with I clusters, consider building the subset randomly. 
First select I nuclei at random for each of the / clusters. 
The remaining k — I nuclei are chosen randomly from 
the neighbors of any previously chosen nuclei. Let L be 
the average number of "nearest neighbors" for each nu- 
cleus. The probability that a nucleus will bridge two clus- 
ters that were meant to be separate is at most O (■^)- 
So the probability that any of the k nuclei will bridge 

two clusters is at most O ( 1 — fl — 



Therefore, as long as k^ ^ N/L, we can accurately es- 
timate the number of subsets in this manner of choosing 
nuclei at random. Under this condition, the number of 
subsets of size k with I clusters will scale roughly (with- 
out dividing out permutations) as 0{N^ L^~^), growing 
exponentially with N as we increase the number of clus- 
ters. Clusters of size one give no contribution; therefore, 
assuming k <^ \/N/L and N :s> L, the possible con- 
tributing subsets of size k are dominated by those con- 
taining all pairs (and a single triplet if k is odd) which 
maximizes the number of contained clusters. If k — 21 
is even, then the number of contributing subsets of size 

k <^ y^N/L, largely composed of I pairs, is O ( ^^-^^ 



\ 2'n 

Therefore, our subset expansion [Eq. H16(l ] error grows 

with fco = 2la < ^/nJL as O ( '^'^i^^^'" ) ■ This is 

problematic because X^NL is not necessarily small. Not- 
ing, however, that most contributing subsets are com- 
posed entirely of pairs (except the one triplet of odd-sized 
subsets), we can approximate the solution up to order 
fco ~ ^/N/L by adding all possible products of pair con- 
tributions as obtained by distributing the product in the 
following, 

VEir) = n {1 + <™W[1 + 0(AL)]} (17) 



-0 



which gives the lowest order of our cluster expansion. 
The correction O ^ ^'^ ^i^C^^ ° ^ is infinitesimal provided 

^0 = fco/2 ^ y/N/L > X^NL, or A < l/iNL^f'^. For 
N ^ 10"' and L ^ 10, we get A < 0.02 as the condition 
for disconnected pairs to dominate spectral diffusion de- 



cay. The 0{XL) in Eq. H17() represents the error incurred 
by not considering clusters larger than pairs (including 
what is required for odd subsets) and can be thought of 
as the contribution you get by adding a neighbor to one 
of the distributed pairs. Not all of the terms obtained 
by distributing Eq. (|17|) will contain pairs that are dis- 
connected from each other. However, using the same 
argument we used to estimate numbers of contributing 
subsets, I < fco/2 random pairs will most likely be dis- 
connected when fco ~ \/N/L. These extraneous terms 
are therefore negligible at each order below fco. 

For a cluster of two nuclear spins (2? = {n, m}) exact 
evaluation of vx>{t) using Eq. (|10l) in the high tempera- 
ture limit leads to 



w«m(r) = l + <™(r) 

= 1 



(1 



[cos {UnmT) - 1] , (18) 
(19) 



with Cnm defined in Eq. H14() . Using Eq. (|17|l and the 
condition yLax{b„m)T *C 1 we write the final expression 
for the lowest order cluster expansion as 



ve{t) « exp • 



V v' 



(r) [1 + 0(AL)]|. 



(20) 



Note that including clusters of three adds a correction 
0{XL) to the decay. 

We have presented two theories. Both require r <C 
Max {bnm)~^\ however, in problems we've considered, the 
decay time is well within this limit. We argued the t ex- 
pansion [Eq. ifT^ ] converges for c„m ^ 1 while the cluster 
expansion [Eq. (|17|) ] for Cnm ^ 1- The cluster expansion 
becomes non-perturbative through the use of the exact 
solution Eq. (|18|l . In fact, in the same way that cluster 
size was related to minimum orders of A in the perturba- 
tion expansion, we can also relate cluster size to minimum 
orders of t in the r-expansion. For example, by taking 
Cnm ^ 1 in Eq. H18|l we recover Eq. (|13|l to lowest order, 
showing that this exact solution interpolates between the 
two regimes at lowest order. For physical problems where 
a wide range of parameters c„m coexist, exact evaluation 
of larger clusters provides a novel systematic approxima- 
tion to the problem of spectral diffusion. 

We use Eqs. (|20|l and (|18|l to perform explicit calcula- 
tions of electron spin echo decay of phosphorus impurities 
in silicon 0, 0, Q ■ In this case the parameter An is given 
by the hyperfine shift of a nuclear spin positioned a vector 
Rn from the donor center. 



87r 



757/ft|*(H„)|' 



(21) 



We used 7s = 1.76 x 10^(s G)-^ and 7/ = 5.31 x 
10^(s G)"^. Here VE'(-R„) is the Kohn-Luttinger wave 
function of a phosphorus donor impurity in silicon, as 
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T [ms] 



FIG. 1: Hahn echo decay ve{t, 9) of a phosphorus donor elec- 
tron spin in silicon, (a) Theory (solid lines) and experiment 
is shown for several orientation angles of the magnetic field 
with respect to the crystal lattice, ranging from the [100] to 
the [110] direction {0 = 0, 10, 20, . . . , 90). (b) Here we plot 
— \nvE{j,6) +\nvE{j,9 = 0) , allowing for the removal of any 
decoherence mechanism which is independent of 9. The quali- 
tative and quantitative agreement between theory and exper- 
iment is remarkable, in contrast to the stochastic approach 
(dashed) . 



described in Ref. 0. The central '^^P nuclear spin does 
not contribute to spectral diffusion because its hyperfine 
energy is significantly larger than any of its neighbors, 
suppressing the spin flips by energy conservation. Dipo- 
lar coupling is given by 



1 2fcl - 3cos^6'„™ 

bn,n - -jllfT' • 



(22) 



It contains an important anisotropy with respect to the 
angle 0nm formed between the applied magnetic field and 
the bond vector linking the two spins (i2„m)- This prop- 
erty leads to a strong dependence of spin echo decay when 
the sample is rotated with respect to the applied B field 
direction. Fig.^shows experimental data when the sam- 
ple is rotated from the [100] to the [110] direction. Here 
the cluster approximation is expected to be appropriate 
for r <C 1 — 5 ms. Finally, in a natural sample of silicon 
only a small fraction / = 4.67% of lattice sites have non- 
zero nuclear spin (these are the spin- 1/2 ^^Si isotopes). 
Averaging Eq. (fT7|l we get 



ve{t)^ n [(l-/') + /V™(r)] 



(23) 



Our numerical calculations of Hahn echo decay for 
several magnetic field orientation angles are shown on 
Fig. n^a). Wc also show experimental data taken for 
bulk natural silicon with phosphorus doping concentra- 
tion equal to 2 x 10^^ cm'^ 8]. The high concentration 



of phosphorus donors leads to an additional decoherence 
channel arising from the direct spin-spin coupling be- 
tween the electron spins that contribute to the echo. This 
contribution can be shown to add a multiplicative factor 
exp (-r/lms) to Eq. ^23 |I2- 

Because this contribution 
is independent of the orientation angle, we can factor it 
out by subtracting the 9 = contribution from the log- 
arithm of the experimental data taken at angle 9. The 
result is shown on Fig. Q^b) (log-log scale). Our theory 
seems to explain the time dependence of the echo quite 
well. To check convergence of our cluster expansion we 
have gone to the next order. Including clusters of three 
amounts to a contribution of 0.1%, in agreement with 
our estimate of XfL ^ 10~^. We have also verified that 
our cluster expansion results agree quantitatively with 
Eq. H13I) for small t when excluding nuclei close to the 
center of the electron wave function where Cnm 3> 1 . This 
result is to be compared with the recent stochastic the- 
ory developed by two of us |5*] (Dashed line in Fig. ^b) 
shows the stochastic calculation for 9 = 60°). Although 
the stochastic theory yields correct order of magnitude 
for the coherence times, it fails qualitatively in explain- 
ing the time dependence. The present method is able to 
incorporate all these features within a fully microscopic 
framework. 

An important issue in the context of quantum infor- 
mation is the behavior of spin coherence at the short- 
est time scales. The experimental data in Fig. ^ reveals 
several oscillating features which are not explained by 
our current method. These are echo modulations aris- 
ing from the anisotropic hyperfine coupling omitted in 
Eq. 1^ 0. This effect can be substantially reduced by 
going to higher magnetic fields (In a quantum computer 
B ^ 9 Tesla will probably be required in order to avoid 
loss of fidelity due to echo modulation (13|j). On the other 
hand, spectral diffusion is essentially independent of mag- 
netic field even to extremely high values {B ^ 10 Tesla). 
Nevertheless this effect is expected to be absent in III-V 
materials hence our theory allows the study of spin 
coherence at time scales of great importance for quantum 
information purposes but currently inaccessible experi- 
mentally. 

In conclusion, we describe a new quantum approach for 
the problem of localized electron spin decoherence due to 
the fluctuation of dipolar coupled nuclear spins. In con- 
trast to former theories, our method requires no ad hoc 
stochastic assumption on the complex dynamics of the 
environment responsible for decoherence. Hence it pro- 
vides an important example where direct integration of 
the environmental equations of motion provides a system- 
atic understanding of the loss of coherence which needs 
to be controlled for quantum information applications. 
We are indebted to A.M. Tyryshkin and S.A. Lyon for 
providing the experimental data shown in Fig. ^ This 
work is supported by ARDA, ARO, and NSA-LPS. 

After the completion of the review process of our 
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manuscript, a preprint by Wang Yao, Ren-Bao Liu, 
and L. J. Sham [l5j appeared exactly reproducing our 
lowest order theoretical result by a completely differ- 
ent technique which treats excitations of nuclear pair- 
correlations as quasi-particles that are non-interacting for 
T <C Max(6„.m)~^- This independent agreement demon- 
strates the validity of our cluster expansion technique. 
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